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We obtain from first principles, i.e. from the quark-gluon dynamics, the Gell'Mann-Ne'eman 
eightfold way baryons in an imaginary-time functional integral formulation of 3 + 1 lattice QCD 
in the strong coupling regime (small hopping parameter k > 0). The model has SU(3)c gauge 
and global SU(3) / flavor symmetries. In the subspace of the quantum mechanical physical Hilbert 
space of vectors with an odd number of quarks, the baryons are associated with isolated dispersion 
curves in the energy-momentum spectrum. The spin 1/2 octet and spin 3/2 decuplet baryons 
have asymptotic mass — 31nfi: and for each baryon there is an antibaryon with identical spectral 
properties. All the masses have the form M = — 31nK — 3k^/4 -|- K®r(K), with r(K) real analytic. 
For each member of the octet r(K) is the same; for each member of the decuplet, r(0) is the same. 
So, there is no mass splitting within the octet, and within the decuplet up to and including 0{k^). 
However, there is an octet-decuplet mass difference of 3k^/4 -|- 0{k'). The baryon and anti-baryon 
spectrum is the only spectrum up to near the meson-baryon threshold of « — 51nK. A decoupling 
of hyperplane method is used to naturally unveil the form of the baryon composite fields (there is 
no a priori guesswork), to show the existence of particles and their multiplicities using a spectral 
representation for the two-baryon correlation. We also obtain the (anti-)baryon dispersion curves 
which admit the representation w{k,p) = —3 In k — 3K^/4-|-ft:"^ Sj=i 2 3 (1 — cosp-')/4-|-r(K, p), where 
r{K,p) is of 0{k^). For the octet, r{K,p) is jointly analytic in k and in each , for small jlmp-'l. 
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Oh; 

^ . In Ref. [ij] a quark model with three flavors {up u, down d and strange s) and an SU(3)/ flavor symmetry was 
introduced to describe hadrons by an eightfold way classification scheme. A dynamical SU(3)c local gauge model of 
quarks and gluons and color dynamics was proposed later, the well-known QCD, as a model for the strong interactions. 
It was shown to be asymptotically free , and perturbation theory was used successfully for high energy phenomena 
K*" but not at low energies. To understand the low- lying energy- momentum (E-M) spectrum and confinement (no isolated 
uarks arc observed) a lattice approximation in an imaginary-time functional integral formulation was introduced in 
The use of this approximation in different contexts, for example the strong coupling expansion, can be found 
e.g. in Refs. [1, [1, Q ■ Numerical simulations on the lattice acquired an important status to determine the particle 
content of the model and to give an answer to other questions which were not attainable using perturbation theory 
(see i). 

^ ^ In a mathematically rigorous treatment, and in an imaginary-time setting, a physical Hilbert space TC and E-M 

\ operators are constructed for the lattice QCD in [1, A Feynman-Kac (F-K) formula is also established. 
■ In a series of papers in Refs. [HI, [H, [11] , we determined the low-lying E-M spectrum of increasingly complex SU(3)c 
lattice QCD models in the strong coupling regime, i.e. with small hopping parameter k > and plaquette coupling 
P = 1/(2(7q) (large glueball mass) obeying /3 ^ k ^ 1; we obtained the onc-hadron and the two-hadron bound-state 
^ ' spectra, up to the two-particle energy threshold. 

Here, we obtain the baryon part of the spectrum in the more realistic SU(3)c lattice QCD model with three flavors, 
in 3 -|- 1 dimensions and in the strong coupling regime. We derive the Gell'Mann-Ne'eman eightfold way baryons 
exclusively from the quark-gluon dynamics. No guesswork is needed regarding the form of the baryon composite 
flelds. We show the existence of 56 baryons and their anti-particles. The baryons have asymptotic mass « — 31nK 
and form the spin 1/2 octet and the spin 3/2 decuplet of the eightfold way. Anti-baryons and baryons have the same 
spectral properties by charge conjugation, and we show that they give rise to all the E-M spectrum in the subspace 
Ti-o C H of vectors with an odd number of quarks, up to near the meson-baryon threshold of « —5 In k. The other 
eightfold way particles are mesons, which lie in the even sector of Ti, and have asymptotic masses « — 21nK. The 
eightfold way mesons are treated in Ref. .14]. 

Our lattice QCD model has the partition function Z = J e'^'^''^'^'^^ dip dtp d^{g) , and for F{'ip, ■0, g), the normalized 
correlations are denoted hy (F) — ^ J F{tlj,ilj, g) e~^^'^'^'^^ dip dtp dii{g) . The gauge-invariant action S = S{ip,'ip,g) 
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is Wilson's action '^l with an SU(3) / flavor symmetry given by 

where, besides the sum over repeated indices a, (3 = 1, 2, 3, 4 (spin), a = 1,2,3 (color) and / = 1, 2, 3 = u, d, s (isospin), 
the first sum runs over u = {u°,u) = {u°,u\u^,u^) & 1^ = {±1/2, ±3/2, ±5/2...} xl?, a = ±1 and ^Ji = 0,1,2,3. 
Here, denotes the time direction and the 3 direction is also called the z-direction. e'^ is the unit lattice vector for 
the ^-direction. At a site u G Z^, tpaafiu) are fcrmionic Grassmann fields (the upper hat meaning the presence or 
absence of a bar) and we refer to a = 1,2 as upper spin indices and a = 3,4 (equivalently, + or — respectively) as 
lower ones. For each nearest neighbor oriented bond < u,u±e^ > there is an SU(3)c matrix U{gu,u±ei^ ) parametrized 
by the gauge group element gu,u±ei^ and satisfying U{guM+ei^)^^ = U{gu+ei^,u)- For simplicity, we sometimes drop 
the U from the notation. To each oriented plaquette p there is a plaquette variable x(f^(5p)) where U{gp) is the 
orientation-ordered product of matrices of SU(3)c, and x is the real part of the trace. M = M{m, k) — m + 2k and, 
given K, TO > is chosen such that M^^ — Sap, so that to = 1 — 2k < 1 in the strong coupling regime. Also, we 

take r**^'' = — /4 ± 7^, where 7*^ = ( ^ ] i 1'' = \ ^ i *n ) i J ~ 1' 2, 3, are the 4x4 Dirac matrices and 



satisfy {7^, j^} — 2(5^^/4, where ^'^'"^ are Pauli matrices. d^{g) is the product measure over non-oriented bonds of 
normalized SU(3)c Haar measures (see Ref. jlSl), and are indistinct integration variables. The Grassmann 
integrals are given in [l^; for k = 0, {tpg-^ (x) ipt^ (jj)) = Sai,a2^aia2Sfif2H^ ~ 2/)j ^'^^ ^^'^ integral of monomials is given 
by Wicks theorem. For more details about the model definition and notation, see Refs. [llill^. For free fermions, 
there is no spectral doubling and the dispersion curve is monotone increasing in momentum p, convex for small \p\ . 

The physical quantum mechanical Hilbert space Ti, and the E-M operators H and , j ~ 1,2,3, are defined as in 
SIM El- Polymer expansion methods [l^, [13] ensure the thermodynamic limit of correlations exists and truncated 
correlations have exponential tree decay. The limiting correlations are lattice translational invariant and extend to 
analytic functions in the global coupling parameters k and f3 = l/(2g§) and also in any finite number of local coupling 
parameters. For gauge-invariant F and G restricted to u'^ = 1/2, we have the F-K formula 

{G,ffffffffF)n = {[Tff^F]QG) , (2) 

where Tq° , Tf , i — 1,2,3, denote translation of the functions of Grassmann and gauge variables hy > 0, x — 
3) G I?, = Tf T| and 9 is an antilinear, order reversing operator which involves time reflection 
[lO|. For simplicity, in Eq. we do not distinguish between Grassmann, gauge variables (rhs) and their associated 
Hilbert space vectors (Ihs) in our notation. As linear operators in 7i, I),=o,i.2,3 are mutually commuting; Tq is self- 
adjoint, with —1 < To < 1, and Tj=i_2,3 are unitary. So, Tj — e^^' defines the self-adjoint momentum operator 
P = {P^,P^, P-^) with spectral points p e = (— tt, tt]'^ and T§ = e^^^ > defines the energy operator H >0. We 
call a point in the E-M spectrum with p — a mass. Also, we let £{X'^, A) be the product of the spectral families of 
To, P\ P2 and p3. 

In order to classify and label the baryon states, we note that the action of Eq. ([T]) has a SU(2) © SU(2) spin 
symmetry at k = in the hopping term, separately in the lower and upper components. Motivated by this property, 
we introduce spin operators Jx, Jy, Jz and = Jx + Jy + Jz acting on the Grassmann field algebra, and which 
obey the usual angular momentum commutation relations. Although we adopt the terminology of the continuum, the 
analogy with the continuum only holds for k = for which the spin symmetry is exact. 

To show the existence of particles up to near the meson-baryon threshold, we obtain a matrix valued two-point 
function G{u, v) which has a spectral representation, derived by using the F-K formula and the spectral representations 
of the E-M operators. The Fourier transform G{p) = J2xei,i- G{x)e~'^P-^ , wherep = {p^,p) are conjugate variables, has 
a lattice Kallen-Lehman type representation which allows us to relate momentum space singularities on the imaginary 
p^ axis to points in the E-M spectrum. We want to show that there are isolated dispersion curves, defining the 
eightfold way baryons and their anti-baryons, up to near the energy threshold of — 51nK. To this end, we consider 
the inverse T{p) = G{p)~^, and show that, for fixed p and k, 

f-i(p) = {cof[f(p)]}Vdetf(p), 

provides a meromorphic extension of G{p) in Thus, the singularities of G{p), or spectral points, are contained in 
the zeroes of detr(p) and the dispersion curves 'w{p) satisfy the equation 

detf(]5° = iw(p),p,K) = 0. (3) 

That T{p)~'^ provides a meromorphic extension of G{p) follows from the faster temporal falloff of T{x = u — v), the 
convolution inverse of G, as compared to G. The faster falloff of T(x) gives us a larger strip of analytic! ty in Tmp'^ 
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which is |Xmp''| < —(5 — e)lnK, < e <C 1, and the analyticity impHes the zeros of detf (p) are isolated, for each p 
and K, leading to the existence of particles in the E-M spectrum. To obtain precisely the umber and behavior of the 
dispersion curves, we need the short distance, low k order behavior of F which in turn follows from that of G. 

To find the appropriate two-point function G, the associated fields that create the low-lying excitation spectrum, 
as well as the falloff of G and F, we use a hyperplane decoupling method. In this method, we consider the correlation 

GlmM=^{L{u)M{v)) , (4) 

for vP ^ , and with L and M containing an odd number of V' fields. For each set of adjacent temporal hyperplanes 
(parametrized by p), which separate the points, we substitute the action hopping parameter n by complex Kp. Con- 
cerning the decay of Glm{u, v), intuitively we pick up a decay factor of Kp for each vanishing Kp derivative at Kp = 0. 
Taking the Kp derivatives at Kp — 0, the r*** derivatives r = 0,1,2,4 are zero by imbalance of fermions and or by 
intcrhyperplane gauge field integration. To analyze the third derivative, we need the gauge integral (see Refs. [6l.[l3j) 
/ 501615026250363^^(5) = £010203 £616263/6- With aU fields at the same point, we let 

^af~ ^o6c "^OQi/i ^"602/2 V'cQ3/3 J 

and use the superscript u {£) to denote that only upper, = 1,2 (lower, = 3, 4) spin components occur. Also, the 
superscript means the coefficient of Kp. For v'^ > vP ^ u° + 1/2 < p < w'^ — 1/2, we obtain 

{L{u)M{v))i^) = E,- s Ji^i^ybifM ^))*°^ Hs(P + 1' ^)Miv))^'^ ~ {LiuW^p, ^))(«) {^.{p + 1, w)M{v))(')] . 

.(5) 

A similar expression holds for m*^ > t;". Note that, with the Levi-Civita's coming from X3, the local, gauge-invariant 
(colorless) fields b^^ and b^jr have naturally made their appearance. Below, we will show that (^"/) basic 
excitations. Their linear combinations create the eightfold way baryons and anti-baryons which are related by charge 
conjugation and have identical spectral properties. 

From the vanishing of the zeroth, the first and the second Kp derivatives, which extends to spatial separations as 
well, we obtain the decay 

|Glm(w, 1^)1 < const (6) 

with |w — w| = — v'^l + lu — v\, \u — v\ = X)i=i 23 — and G{p) is analytic in the strip |2rmp°| < — (3 — e) ln«;. 

We emphasize that there is no guesswork regarding the form of the baryon fields in Eq. ([S]), and in the sequel we 
only consider these fields. The appropriate choice for the two-point correlation in our method is to choose L and M 
so that the correlation on the Ihs of Eq. ^ is the same as those on the rhs, i.e. to have closure. For this, we take 
L — 6 and M = b^j^ with lower indices in a, such that only the first term in Eq. ([5]) survives (choosing L — b"j. , 

M = &~ only the second term survives and leads to antibaryons) . Now, with a view to obtaining a convolution 

inverse with a faster temporal decay, we note that there are redundancies in fe^^, which prevents us from defining 
the inverse of Glm- The linear dependencies are eliminated using what we call the totally symmetric property (tsp) 
which is invariance of b under the exchanges fi ^ aj fj . By considering only one element per equivalence class we 
have only 56 independent fields, instead of the original 6'^ = 216 (see Ref. [18]). The number 56 arises as it is the 
dimension of the totally symmetric subspace of the 3— fold product of a 6-dimensional vector space. The dimension 
6 comes from 3 (flavors) x 2 (lower spins). With this restriction in effect and without changing notation we can now 
introduce the normalized fields B^^= 6^y/[ngy] such that, for coincident points, {B^j-B^, p)'^^^ = —Sss'^fji- To see 
the last equality we use the basic formula [b^j-b^, pY^^ = —6 perm(^), where perm is the permanent (similar to det, 
but with only -|- signs), and Aij = (5q,.q' ^z^/', for i,j = 1,2,3. 

Taking into account the other time ordering, we define the two-baryon function for all u and v by (x is the 
characteristic function here) 

Gi^i^{u,v) ^ {Bi^{u)Be^{v)) XuO<vO - {Bi^{u)Bi^{v))* Xu«>vO , (7) 

where now the £'s are collective indices for cJ/'s, and we suppress the lower spin superscripts. For u'^ v^, 

^^2 " ~ X! ^^3^""' (P''^))G'£3L((P+ l''^)'«)X«o<i;0 - X! ^fX^"^' (P + l,w)) GfJf,^{{p,w),v)Xu'>>vO , 



which we write symbolically as 

G(3)(u,w) = -[G(") oG(°)](m,v), 



(8) 
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with true convolution in space. This important relation, which we call the product structure, is instrumental in 
showing the faster temporal F decay, as it feeds into the formula for the 3rd Kp derivative of F at k = 0. Indeed, 
letting G = Gd + Gn, where Gd^t^t^ {u, v) = Ge^t^ (u, u)Se^e2^uv is the diagonal part of G, we define F by the Neumann 
F = (Grf+G„)-i = J2T=o Gd^i-GnG^^f, which converges by the global boimd on G of Eq. ©. Using FG = 1 = GF, 



the Leibniz formula d'^V = X]I=o 



T'd'-'GdT' (with F' = -F), the product structure of Eq. dH) and the Kp 



derivatives of G, we obtain F(''"°'^'^^(m, v) = 0, for 
and, by imbalance of fermions, T^^\u,v) — 0, - 
follows the faster decay bound 



.[r(o)G(3)F(o)](u,v) = 0, if > 1 

v^\ > 1. From this, using joint analyticity and Cauchy estimates. 



u'^-v^l > 1, F(3)(m,w) 



\Te,e,{u,v)\ < const \k\^ |^|5(|n°-.°|-i)+3|a-iT| ^ |^ 



> 1, 



(9) 



and the rhs is replaced by const if = . Hence, f (p) is analytic in the larger strip |Xmp°| < — (5 — e) Iuk. 

To relate points in the E-M spectrum to singularities of G i-^^ir^ijp) , we first use the F-K formula to obtain a spectral 
representation, with Bi = Be{l/2, 0) and x = v — u, 

GiAi^) = -{Bi„f\'-"\f^Bi,)n = ~ [ I {X° f^-^e-'^-^dxiBi„£{\°,X)B,,)n, 

for X e Z'*, x° ^ 0, and is an even function of x by parity symmetry. For the fourier transform, after separating out 
the a;*^ = contribution, we get 



(10) 



with f{x,y) EE (e" - y)"^ + (e"" - where dAoap.^i^sC^") = /t3 ^{p- ^) d),odx{Bij^,£{X° , X)Bi^)n, and we have 

set G{p} = J2x e~'P-^G{x° ^0,x). 

Singularities on the Tm p^ axis are spectral points and are contained in the zeroes of det T{p). We first restrict our 
attention to the determination of the masses (i.e. p = 0), which is simplified passing to a basis where T{p'^,p = 0) is 
diagonal. The diagonalization is achieved by fully exploiting the SU(3) / symmetry, and passing to the eightfold way 
baryon particle basis. The particle basis is related to the individual spin and isospin basis we have dealt with up to 
now by a linear real orthonormal transformation. This transformation preserves the product structure, and thus the 
larger strip analyticity region of Tj.-^r2J where now the r^'s are collective indices of the particle basis. 

The use of the flavor symmetry reduces F(p) to a block form with 8 identical 2x2 blocks associated with the spin 
1/2 octet, and 10 identical 4x4 blocks associated with the spin 3/2 decuplet. The octet basis is given by 



P± - 
■=■0 



3V2 



^{-^a+u^b-d - '^a+d'^b-u)'^c±d, 



{■tpa+u^b- 
S± = "^{^a+u^b- 



3V2 ^ 



Uc±s , =± = f^i-^a+dlpb- 



ila+si>b-d)i', 



cibs 7 



and the decuplet basis is given by 



V'a+s'06-«)V'c 

'lpa+s^b-u)^l^c±u , T,'^ ^ ^i2il:a±ui>b±di>cTs - 4'a-u1pb+d^c±s 
4^a+s'ipb-d)i'c±d , ^± ^ ^{'ipa+u'ipb-d - 'il'a+d'4'b-u)-ii^c±s, 



1pa+ui^b-di^c±s), 



Ati = 

2 

2 

Ail = 

2 

A++ - 

v*+ 
^ ±3 



^ ±3 
S ±3 



^i^a±u^b±u^cTd 



2lpa±u'ipbTui'c±d) 
^{'^■'Pa±uA±d1pcTd + i^a^ui'b±di^c±d) 
1^ 4'a±d4>b±di'c^d 



2^3 



I— >* — 

^ ±3 

2 

+3 = 



±3 ~ 6 



^^i^a±u1pb±u'4'c±8 
^ 'ipa±uA±d^c±s 
^4'a±dTpb±d'4>c±s 
^tpa±ui^b±s'tpc±s 
^'ipa±di'b±s'ipc±s 
^^a±s'ipb±s'4>c±s 



V*0 

2 

" ±1 
2 

^ ±1 



■lpa±uA±ui^c±d, 
^■fpa±u-ipb±d-ipc±d, 
^ i>a±d1pb±d1pc±d, 
^ 1pa±u1pb±ui>c±u, 



A|3 

2 

AI3 
Aft 



■^{lpa±u'ipb±di^cTs + i^a±utpbTdi^c±s + i^aTutpb±d1pc±s) , 



Cob 

6 

lob 

6 



f^{-4'a±d'ipb±d1pcT^^ + '2->Pa±d->PbTd'ipc±s) 
{i^a^di^b±s + '2'ipa±d'ipbzf:s)i>c±s, 
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We use the same particle symbol on the left, as for the Gell-Mann-Neeman eightfold way in the continuum, and the 
barred fields that creates it at the right. The particle basis collective indices specify the quantum numbers of total 
isospin /, third component I3, the total hypercharge Y and the value of the quadratic Casimir C2 {C2 = 3 for the 
octet and 6 for the decuplet.), besides the labels of total spin J and the z-component J^. Jz is theparticle symbol 
subscript. These labels are the same as the conventional ones for the continuum found in e.g. Refs. P, U^. For fixed 
Jz, the octet (decuplet) vectors form a basis for the 8 (10) dimensional irreducible representation of SU(3)/. Here 
we naturally obtain one 8 and one 10 dimensional representation of SU(3)/, which is not an ad hoc choice as in the 
3x3x3 = 10 ©8©8©1 decomposition appearing in the non-dynamical group theoretical construction. 

The j-th component of total isospin, for j — 1, 2, 3, is defined through the linear operator acting on a function F of 
the Grassmann algebra by AjF = lime\o{[-F({t'iV'}, {i^U^j}) - FiW, where Uj = Uj{e) = exp(iAj6'/2), 

j — 1, ...,8 is an element of SU(3)/ and the \j are the usual Gell-Mann matrices P, [l^. Y is defined as 2As/\/^ and 
C2 ~ 8^j- Defining the linear operator W{U)F — i^({[/-0}, {V'f^^}) then W{U) lifts to a unitary operator 

W{U) on H. by using the F-K formula and the SU(3)/ symmetry. The generators Aj = \img\^o{W{Uj) — l]/{i9)} 
of the eight one-parameter subgroups are self-adjoint operators in H. I3, P = -I- /| -I- /|, F, C2 are mutually 
commuting and their eigenvalues are quantum numbers which are used to label the states. Jx, Jy, and Jz are also 
defined similarly on the Grassmann algebra only, with [/ = C/2 ® U2 = exp{iaW /2) G SU(2). Of course, formally 
taking the continuum limit, the generator Jz is the usual z-component generator for rotations of spin. The eigenvalues 
of Jz and are also used to label the states. 

We now use the auxiliary function method to determine: the octet mass (they are all the same), the decuplet 
masses (they are all the same up to and including 0{k^)), and the 0{k^) octct-decuplet mass splitting. By the use 
of the symmetry of 7r/2 rotations about e^, the matrix T{p) is seen to be diagonal at p = 0; and by reflections the 
elements only depend on \Jz\- The determinant factorizes, and we consider one of the 56 typical factors (for which 
we omit all indices). As the mass runs out to infinity as k \ 0, the usual implicit function theorem does not apply 
to solve Eq. ([3]) at p = 0. We make a nonlinear transformation from p'^ to an auxiliary variable 

w -1~ C3ip)K^ + K,^e-'P\ (11) 

with C3{p) = — 2 3 cosp^/A, and introduce an auxiliary function H(w,k) such that r(p°,p) = H[w = —1 — 

C3{p)k^ + K'^e^'^", k). By doing this, we bring the non-singular part of the mass M -I- 31n«; from infinity to close to 
w = 0, as K \ 0. Using time reversal and parity symmetry, we have r(a;°, x) = r(— x", x), and i?(w, n) is defined by 

i7(w,K) = ^r(a;° = 0,f)e-'P-^+ ^ T{n,x) 

X x,n—1.2... 

The global bound on F of Eq. ^ guarantees that H{w, n) is jointly analytic in k and w. However, we want to control 
the mass to order and for this we need the low n order short distance behavior of r(a;), which in turn follows 
from that of G{x). More precisely, we need F(a;'^ = n,x) j tc'^) up to and including order k®. The normalization 
condition G^'^\x = 0) = —1 implies F'-'^)(a; = 0) = —1 and, by a simple argument, the product formula gives 
G{x — e^) — -f- 0{k^), which implies F(x = e'') = k'^ + ©(k"*). Other contributions are found by explicit 
calculation of coefficients of the hopping parameter expansion of G{x). Namely, there are contributions arising from 
non-intersecting paths connecting the point to 2: and paths that emit and absorb a meson. After a lengthy calculation, 
we find H{w, k) = w + + uqk^ + bn^ + X]n=i 4 '^'sn+ei^ + + h{w, k)k'^ , with b and the c's taking the 
same value for the octet and the decuplet, and h(w,K) jointly analytic in w and k. Note that, in the term aeK^, we 
have separated the contributions coming from all points of the form x = ee^ + e' , with ij = 12, 13, 23, e, e' = ±1, 
which we call spatial angles, ag takes the value Co ~ 3/8 (a^ = ~3/8) for the octet (decuplet). As 7?(0,0) = and 
[9i//9w](0, 0) = 1, the analytic implicit function theorem implies that H{w,k) ~ has the analytic solution w{k) 
given by, with b' = b+l + E„=i,...,4 ^3,1+6' 

w{k) ^-ae-b'n^ + 0{k'^). 

Returning to Eq. (fTT|) and setting p'~' — iM, we get 

M = - In + + C3{6)k^ + w{k,)) = -31nK + C3(0)K^ - [ae + b' + 03(0)^ /2]k^ + 0{k'^) , 

and the octet-decuplet mass difference is Md — Mo — (oo — adjn^ + C>{k^) = 3k^/4 + ©(k^). 

For the Py^O dispersion curves, the 2x2 and the 4x4 blocks of T{p) still have a complicated structure even after 
the use of the usual well known symmetries. However, we have found a new local symmetry of spin flip J-s which 



C3{P)K 



C3(p)K^ 



-ip.x 



(12) 
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is a composition of the nonlocal, antilinear time reversal, the local, linear charge conjugation and a nonlocal, linear 
time coordinate reflection ^] . This symmetry shows that the 2x2 octet blocks are diagonal and a multiple of the 
identity. For the octet, the identical dispersion curves w{p) can be obtained by using the auxiliary function method 
as done before for the masses. For the 4x4 decuplet blocks, the spin flip symmetry simplifies the matrix T{p) but 
it is not diagonal. We have not been able to apply the auxiliary function method. However, we can use a Rouche 
theorem argument (principle of the argument) on detT4x4{p) to show that for each fixed p there are exactly four 
pairwise identical solutions. 

The dispersion relations (curves for the octet) admit the representation w{k,p) = — 3 In k — 3k"^/4 + X]j=i 23(1^ 
cosp^)/4 + r{K,p), where r{K,p) is of 0{k^). For the octet, r{K,p) is analytic in k and in each , for small |Xmp^ |. 

Up to now we have shown the existence of baryons in the subspace Tib C Tio, the subspace generated by vectors of 
the form b^j: The E-M spectrum up to near the meson-baryon threshold of — 51nK consists of dispersion relations 
associated with the 56 eightfold way baryons and the 56 anti-baryons. In principle, there may be other states in Ho 
with spectrum in the same energy interval. We show this is not the case using a correlation subtraction method Tlj . 

We make some concluding remarks. We have considered a model with the same quark mass and flavor SU(3)/ 
symmetry, which has baryons and antibaryons. Our method extends to treat rigorously the case of SU(7V)/ to 
uncover {2N + 2)!/[3!(2A^ — 1)!] baryons and also to treat the case with different quark masses and broken flavor 
symmetry. The product structure still holds, the elementary excitations are revealed and are the same as the b^^ 
obtained here. The difficulty in obtaining the dispersion curves for the decuplet particles disappears for the continuum 
if the decuplet fields transform under the Poincare group. In this case, the four identical dispersion curves are of 
course the relativistic ones. 

The determination of the one-hadron spectrum is a necessary step to analyze the existence of bound states. Our 
method is powerful enough to access the hadron-hadron spectrum and should help in clarifying fundamental open 
questions as e.g. the existence of certain meson-baryon bound-states and pentaquarks. 

This work was supported by CNPq and FAPESP. 
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